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The flavor-changing neutral current process b — )■ s££ 
occurs in the standard model (SM) through higher order 
diagrams where new physics contributions could arise. 
Accurate SM predictions make the b -> s££ phenomenol- 
ogy suited to uncover early indications of new physics 0- 
[1] , especially through observables like the lepton forward- 
backward asymmetry (Afb) and the differential branch- 
ing fraction (BR) as a function of dilepton mass Mu. The 
6 — > s££ amplitudes can be described in terms of short dis- 
tance operators and effective Wilson coefficients Cy^g^io- 
Some new physics models [l[ allow the flipped sign of Cy. 
This results in the opposite sign of Afb in the small 
region (g2 = M|^c2). Recently, BaBar and Belle ^ mea- 
sured an ApB in the i?" — K*^£'^£~ decay larger than 
the SM expectation. The decay S° (j){1020) fx+ fx' has 
not been seen in previous searches by CDF Q and DO Q . 

In this Letter we report an update of our previous 
analysis Q of the rare decay modes K^fi^ii^, 
— )• K*°iJ,'^fi~, and — 4'fi'^fi~ using an increased 
data sample of pp collisions at a center-of-mass energy 
of y/s = 1.96 TeV corresponding to an integrated lumi- 
nosity of 4.4 fb~^, collected with the CDF II detector 
between March 2002 and January 2009. We update the 
BR measurements and also report the measurement of 
ApB in the B° K*°fi'^fi~ decay. 

We reconstruct B — >■ hfjC^ji'^ candidates, where B 
stands for B+, B°, or and h stands for K+, K*°, or 
(j), respectively. Charge-conjugation is implied through- 
out the Letter. The K*'^ {(j)) meson is reconstructed in 
the decay — >■ K+tt" {(j) — > K+K^). We also recon- 
struct B — )• Jf-tph decays as normalization channels in 
BR measurements, because they have final states identi- 
cal to those of the signals, resulting in a cancellation of 
many systematic uncertainties. The relative BR's are: 



BR{B J/iJ)h) 



N 



cBR(J/V' ^ 

(1) 

where Nh^,+^- (Nj/^h) is the B h^i+ n {B J/iph) 
yield, and E/i^+^i- /£j/^/i is the relative reconstruction ef- 
ficiency determined from the simulation. 

The CDF II detector is described in detail in Ref. fzl 
with the detector subsystems relevant for this analysis 
discussed in Ref. Q. 

A sample of dimuon events is selected by the on- 
line trigger system. The trigger requires two opposite 
charged particles with transverse momentum > 1.5 or 
2.0GeV/c depending on the trigger condition, matched 
to the muon chambers. We use the muon chambers de- 
tect muons within \r]\ < 0.6 and 0.6 < It;] < 1.0 "9^. 



The trigger also requires L^y > 200 /xm, where L^y is 
the transverse displacement of their intersection from the 
beamline. The detail of the trigger system and selection 
requirements can be found in Ref. [5||. 

The offline loose event selection begins by looking 
for a common vertex of two trigger muons with one 
(two opposite-charge) reconstructed charged particle(s) 
to form a B+ ^ K+ n+ ^i,- (a B° -J> K*°^jL+^i- or a 

The probability of the ver- 



) candidate. 
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tex fit is required to be greater than 10 All 
charged particle trajectories are required to be associ- 
ated with hits in the silicon vertex detector and to 
have pt > 0.4GeV/c. In addition, we require 
PT{h) > l.OGeV/c and priB) > 4.0GeV/c. We require 
that the B candidate's decay is consistent with being dis- 
placed from the primary interaction point in the trans- 
verse plane by Lxy{B) /a{Lxy{B)) > 3, and a{Lxy{B)) 
is the estimated uncertainty of Lxy{B). We also require 
that the B candidate comes from the primary vertex by 
Mo(^)| < 120 ^m, where do{B) is the distance of closest 
approach of the B trajectory to the beamline. 

For B° (B°) candidates the K+n- {K+R-) mass 
must lie within 50 (10) MeV/c^ of the world average 
(0) mass [13] • The ambiguity of the mass assignment in 
the — ^ if +7r~ decay is handled by choosing the com- 
bination with the if +7r~ mass closer to the known if *° 
mass. This results in the correct mass assignments for 
about 92% of the decays as determined from the simula- 
tion. Particle identification is performed with the timc- 
of-ffight and the ionization energy loss (dE/dx) probabil- 
ities of the particle hypothesis. We require loose parti- 
cle identification for both kaons and pions coming from 
the if *° meson or (j> meson to reduce combinatorial back- 
ground. This removes 15% of the B mass sideband events 
while 99.5% of the signal is retained. We also require a 
muon likelihood ll| to suppress hadron tracks that pro- 
duce false trigger muons. 

Rare decay candidates with a dimuon mass near 
the J/ip {ip') are rejected: 8.68 (12.86) < q'^ < 
10.09 (14.18) GeVVc^- To eliminate the radiative 
charmonium decays that escaped rejection above, we 
remove candidates consistent with originating from 
a, B ^ J/xjj{')h decay followed by the decay 
of the J/rp{') into two muons and a photon: 
\iMi^^^,h)~M^^^)~iMi^^^^)-Mff^^^,^)\ < 100MeV/c2, 
where the PDG superscript indicates known experimen- 
tal averages im] and M{fifj.) < M^f^^,y We also reject 
candidates if an opposite sign hadron-muon combination 
of the daughters, assigned the muon mass, satisfy J/^p 
or t/y mass within 40MeV/c^. This removes charmo- 
nium decays where one of the muons is misidentified as a 
hadron. We reject candidates in which two (three) track 
combinations are compatible within ±25MeV/c^ with 
D° K-TT+ {D+ if-7r+7r+ or D+ K+K-TT+) 
decays for i3+, B^ ^ and B^ decays, respectively. This 
removes B — > Dtt {D — , and D+) decays where 
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two hadrons are misidentified as muons. 

We train an Artificial Neural Network (NN) classifier 
on simulated signal and a sample of events representative 
of the background events in the signal region. To simu- 
late the signal we use PYTHIA and EvtGen [l^ based 
on the SM expectation [ij. The background sample is 
obtained from the sidebands of the B invariant mass dis- 
tribution. We take only the higher mass sideband for the 
and B° decays since the lower sideband is populated 
with physics background from partially reconstructed B 
meson decays. We use both sidebands for B^ decays. We 
use 7-10 observables based on B and daughter's kine- 
matics (e.g. pt and mass), vertex qualities, and muon 
likelihoods. We optimize the NN threshold in order to 
maximize both the BR and the Afb significance. For 
the B^ and B° analysis we optimize the NN threshold 
by maximizing Ns/\/Ns + iVf,, where Ns (Nt,) is the ex- 
pected number of signal (background) events. We deter- 
mine Ns by Eq. ([T]) with the world average BR and NN 
cut efficiency of the simulated signal events, and deter- 
mine Nb from the number of sideband events scaled to 
the signal region, which is defined as ±2a from the world 
average B mass, and NN cut efficiency of the sideband 
events. For decays, Ns is taken from a theoretical 
prediction [TT 
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We maximize Ns/{5/2 + y/Nt) 
The signal yield is obtained by an unbinned max- 
imum log-likelihood fit to the B candidate invari- 
ant mass distribution. The likelihood is given by 
^ = n(/sig'^sig + (1 - /sig)^bg), where is the sig- 
nal fraction, Vsig is the signal probability density func- 
tions (PDF) parametrized with two Gaussian distribu- 
tions with different means, and Vhg is the background 
PDF modeled with a first- or second-order polynomial. 
The signal PDF's are determined from the simulated sig- 
nal and the B mass resolution is scaled by the ratio of 
the mass resolution in J/tph data and simulation, which 
ranges from 1.07 to 1.09. The background PDF's are 
determined from sideband data. Fitted parameters are 
/sig, the mean B mass, and the background shape. The 
fit range for B+ and B° {B°) decays is from 5.18 (5.00) 
to 5.70GeV/c^, to avoid the region dominated by the 
physics background. 

While the contribution from charmless B decays 
is negligible due to the muon identification, we find 
a sizeable crosstalk between K*'^^^ii^ and 

B^ — ^ (j))!^ IJL~ contributing approximately 1% of the sig- 
nal, as estimated from simulation. These contributions, 
whose fractions are determined by simulation assuming 
the world average BR and the theoretical prediction [l^ , 
are subtracted from the fit results for the signal yields. 

By optimized NN threshold we reject 99.5-99.8% of 
background events in the signal region. Figure [1] shows 
the B mass distributions. The statistical significance is 
s = \/ —2 ln(£nuii/'Cmax), whcrc £max is obtained from a 
fit with the signal fraction free to fioat and the mean B 
meson mass fixed to the fitted value in the corresponding 



o :Yield:101±12 

^ 607Mass:5284 ± 3 MeV/c^ 
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""b I'.l — Signal 
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- Data 
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FIG. 1: Mass of S" -> A'*V+M" and B° -> can- 
didates with fit results overlaid. The vertical lines show the 
signal region. 



normalization channel, and >Cnun is the maximum likeli- 
hood obtained from a fit with /gig = 0. Systematic un- 
certainty is not considered in the significance evaluation. 
We obtain s = 8.5a, 9.7cr, and 6.3cr for B+, B", and 
decays, respectively. The observed yields are listed in 
Table m This is the first observation of the — (j)fi^fj,~ 
mode. 



TABLE I: Summary of observed yields. The numbers in 
parentheses are the number of events in the signal region. 



Mode 



N 



J/jih 



IT 121 ± 16 (218) 43704 ± 245 (55296) 0.434 ± 0.006 
B° 101 ± 12 (140) 15815 ± 178 (22952) 0.477 ± 0.009 
Bl 27 ± 6 (40) 2930 ± 64 (3883) 0.498 ± 0.012 



We do not apply a NN selection to J j-^h channels, be- 
cause these signals are of sufficient size and purity with 
the loose selection. To obtain the relative efficiency of 
Eq. ([1]), the NN cut efficiency of the loosely selected 
events is considered in addition to the relative efficiency 
of the loose selection. 

The dominant source of systematic uncertainty for 
each BR measurement is the background PDF param- 
eterization (3.9%) for i?"*", the discrepancy of the NN cut 
efficiency between data and simulation (4.8%) for 
and particle identification (3.5%) for B^. For the abso- 
lute BR measurements we assign the uncertainties of the 
world average BR(B J/i^h) [13]. 

Results of the relative BR (Eq. ([1])) measurements are 
listed in Table |lll The BR statistical uncertainties in- 
clude the Poisson term from finite statistics of the sam- 
ple. We also show the absolute BR which is obtained by 
replacing the normalization channel's BR with the corre- 
sponding world average ^0|] value. 

These numbers are consistent with our previous re- 
sults [1] , B-factory measurements 0, [l^ , and theoretical 
expectations 13| . We also measure differential BRs with 
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TABLE II: Measured branching fractions of rare modes. First 
(second) uncertainty is statistical (systematic). 



Mode 



Relative BR(IO^) Absolute BR(10' 



B+ K+fj,-* 



fi- 0.38 ± 0.05 ± 0.02 
II- 0.80 ±0.10 ±0.06 
1.11 ±0.25 ±0.09 



0.38 ±0.05 ±0.03 
1.06 ±0.14 ±0.09 
1.44 ±0.33 ±0.46 



respect to the dimuon mass. Events in the signal mass 
region are grouped into independent bins. Figure [2] 
(a, b) shows the differential BR for i?+ K^fi^fj,- and 
B° K*°n+ii-. 

The ApB and the K*'^ longitudinal polarization frac- 
tion (Fl) are extracted by an unbinned likelihood fit to 
the cos9fi and cos Ok distributions, respectively, where 
6^ is the angle between the /i+ (/J^) direction and the 
direction opposite to the B (B) meson in the dimuon 
restframe, and 9k is the angle between the kaon direc- 
tion and the direction opposite to the B meson in the if *° 
rest frame. The differential decay rates [3] are sensitive 
to cos 9k and cos 9^ through the angular distributions 
given by |Fl cos^ 6'k + |(1 - Fl)(1 - cos^ 9k) for cos 9k 



and|FL(l- 



- cos^ ) ± I ( 1 - Fl ) ( 1 ± cos2 6*^ ) + ApB cos 61^ 



for cos0^. We measure Fl and Ape for B° K*°fj,'^fi~ 
and also Ape for B^ — > K^fi^fi^. Angular acceptances 
are obtained from simulated signal samples assuming un- 
polarized decays. 

The contribution from decays with K-n swapped K**^ 
mesons distorts the signal distribution and swaps the sign 
of cos 6*^. This effect is considered by adding an addi- 
tional signal-like term to the likelihood function. The 
contribution from decays with non-resonant K-ir is con- 
sidered to be small [2] and neglected in the fit. For the 
5+ decay, we set Fl = 1 and consider no scalar term Q . 

The combinatorial background PDF shape is taken 
from the B mass upper sideband that is used for the 
NN training. In the fit to cos Ok (cos 6*^) distribution, 
the only free parameter is Fl (Afb)- For the cos 6*^ fit, 
the value of Fl is fixed to the cos 9k fit result. 

Most dominant source of systematic uncertainty for 
each angular fit is the fit bias near the physical bound- 
ary (0.02-0.07) for Fl in the uncertainty of the Fl fit 
(0.02-0.12) for Afb in B^, and the angular background 
shape (0.01-0.07) for Afb in B'^ . The total systematic 
uncertainties lie in the range 0.02-0.08 for Fl in 0.05- 
0.25 for Afb in S", and 0.02-0.08 for Afb in B+ . The an- 
gular fit results are shown in Fig.[2](c, d) and summarized 
in Table Uni Results in the range < 9^ < 4.3GeVVc^ 
and 1 < g2 < GGeVVc^ are also included. 

In summary, we have updated our previous analysis 
of the flavor-changing neutral current decays b — >■ sfifi 
using data corresponding to an integrated luminosity 
of 4.4 fb~^. We report the first observation of the 
— >■ the most rare B^ decay observed to 
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FIG. 2: Differential BR of B+ K+ fi+ (a) and dif- 
ferential BR (b), longitudinal A'*" polarization fraction (c), 
and forward-backward asymmetry (d) of B° — >■ K*^^ fi'^ fj,~ , 
as a function of squared dimuon mass. Points are the fit re- 
sult. The solid curves are the SM expectation [l|]. Two solid 
curves in (a, b) use maximum- and minimum- allowed form 
factors on differential BR plots. The dotted curves are the 
Ct — — Cf^ expectation. The dashed line is the averaged 
expectation in each squared dimuon mass bin and hatched 
regions are charmonium veto regions. 



date, and measure the total BR. We measure the to- 
tal BR, differential BR, Afb of the B+ K+ ^1+ ^~ and 
B^ K*'-* jj.^ fi^ , with respect to q^. We also measure Fl 
of 5° — K*'^fi^fi~ prior to Afb- These are consistent 
and competitive with the other current best results. At 
present there is no evidence of discrepancy from the SM 
prediction. 
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TABLE III: Summary of B° K*° (1+ fT and B+ K+n+fj,' fit results. Maximum q'^ is 19.30 (23.00) GeV^/c^ for B° 



BR(g") (10^ 



Fl(B") 



(B^ 



BR(g+) (10'' 



(GeVyc-') 



-053^ 

0.82t»:- 

,+0.19 
-0.18 
U.OO_Q^8 

0.091^:1! 



ApB 

-0.75 
f0.40 
-0.41 

U.UU_o 28 

O+0.23 
-0.20 



Afb(B+) 



-0.15^ 



II ± 0.08 
± 0.07 
28 ±0.03 
Is ±0.07 
11 ± 0.04 



[0.00, 2.00) 

[2.00, 4.30) 
[4.30, 8.68) 
[10.09, 12.86) 
[14.18, 16.00) 
[16.00, 19.30(23.00)) 



0.98 ± 0. 
1.00 ±0, 
1.69 ±0. 
1.97 ±0. 
1.51 ±0. 
1.35 ±0. 



40 ± 0.09 

38 ±0.09 
57 ±0.15 
47 ±0.17 

36 ±0.13 

37 ±0.12 



0.31] 



0.07 

0.08 
0.07 
0.02 
0.02 
0.03 



0.13 
0.19 



7+0.23 



0.66^ 
0.42l»:« 



0.25 
0.14 
0.05 
0.07 
0.09 
0.10 



0.38 ± 0. 
0.58 ±0 
0.93 ±0 
0.72 ± 
0.38 ± 0. 
0.35 ± 0. 



16 ± 0.03 

19 ±0.04 
25 ± 0.06 

17 ±0.05 

12 ± 0.03 

13 ± 0.02 



+U.31 
0.33 



+0. 



+0. 

0, 
+0. 

0, 

0.03l^ 
0.071°; 



0.72 
-0.20 
-0.10 



TUT! 

0.; 



[0.00, 4.30) 



;± 0.55 ±0.18 



0.03 



0.21 



±0.05 



0.96 ±0.25 ±0.06 



0.36 



±0.06 



+U.2Y 



+U.3b 



0.: 



[1.00,6.00) 



1.60 ±0.54 ±0.14 0.50 



±0.03 



0.43 



±0.06 



1.01 ±0.26 ±0.07 



0.08 



±0.07 
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